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ABSTRACT 

We present a model for the rotational evolution of a young, solar mass star interacting with an 
accretion disk. The model incorporates a description of the angular momentum transfer between the 
star and disk due to a magnetic connection, and includes changes in the star's mass and radius and 
a decreasing accretion rate. The model also includes, for the first time in a spin evolution model, the 
opening of the stellar magnetic field lines, as expected to arise from twisting via star-disk differential 
rotation. In order to isolate the effect that this has on the star-disk interaction torques, we neglect the 
influence of torques that may arise from open field regions connected to the star or disk. For a range 
of magnetic field strengths, accretion rates, and initial spin rates, we compute the stellar spin rates 
of pre-main-sequence stars as they evolve on the Hayashi track to an age of 3 Myr. How much the 
field opening affects the spin depends on the strength of the coupling of the magnetic field to the disk. 
For the relatively strong coupling (i.e., high magnetic Reynolds number) expected in real systems, 
all models predict spin periods of less than ^ 3 days, in the age range of 1-3 Myr. Furthermore, 
these systems typically do not reach an equilibrium spin rate within 3 Myr, so that the spin at any 
given time depends upon the choice of initial spin rate. This corroborates earlier suggestions that, 
in order to explain the full range of observed rotation periods of approximately 1-10 days, additional 
processes, such as the angular momentum loss from powerful stellar winds, are necessary. 
Subject headings: 97.10.Cv, 97.10.Gz, 97.10.Kc, 97.21. -Ha 
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1. INTRODUCTION 

The last three decades of observational surveys of pre- 
main-sequence stars have resulted in the measurement 
of rotation periods and/or rotational velocities (w sin^ 
of thousands of stars (see reviews by Herbst et al. 2003; 
lScholzl[2009[ ). The analysis of these data have revealed 
a number of mysteries regarding the spin rates of young 
stars. Perhaps the most promin ent of these mysterie s 
is that which was first noted by iVogel fc Kuhil (|1981[ ). 
Specifically, a large fraction of nearly solar mass pre- 
main-sequence stars rotate much more slowly than ex- 
pected. The fact that stars are built up from material 
with high specific angular momentum, and that these 
young stars are still contracting, leads to the expecta- 
tion of a rotation rate near the breakup velocity. 

Solar mass stars with ages less than a few Myr have 
rotation periods typically in the range of 1-10 days but 
with a tai l in the distribution extending to around 20 
davs fe.g..lAttridge &: HerbslJll992HChoi fc Herbsti | 1996 : 



Stassun et al.l 119991: iRebuUT [20011: i RebuU et al.n 2004: 
Covev et al.ll2005HScholdl2009[ ). The statistics show that 
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about half of these stars are indeed fairly rapid rotators 
and do seem to spin up as t hey approach the main sc- 
guence (iBouyier et al.l 119971: Ide la Reza fcTInzon 2004; 
iScholz et all 120071: llrwin et al.ll2008[) . However, roughlv 
half of the stars younger than a f ew Myr rota.te at 
about 10% or less of breakup sp eed (jRebull et al.ll2003 : 
iHerbst et all 120071 : iScholzl l2009( ) . Thus, there is some 
mechanism operating that is capable of removing sig- 
nificant amounts of angular momentum during the pre- 



main-sequence phase. 

Following the s uggestion oflEdwards et al.l (I1993D, the 
modehng work of lBouvier et al.l ()1997[ ) and lRebuU et al.l 
(|2004il showed that the observed range of spins could be 
reproduced reasonably well, by assuming that the pres- 
ence of an accretion disk somehow results in a constant 
stellar spin period of around one week. There is some 
observational evidence that the population of stars with 
disks on average rotat e more slowly than those with- 
out disks (e.g. , [Edwards et al.lll993t[Choi fc Herbstll 19961: 
Stassun et all ri999; .Herbst et al.l I2OOOI: iLittlefair et al.l 






20051: IRebull et al.l boOfit ICieza fc Baliberi 120071 ). Also 
the idea that the angular spin rate is held at a constant 
by the presence of the disk is loosly based on physical 
models for the interaction between a star and surround- 
ing disk. 

There are generally two prominent theoretical ideas 
for how accreting stars can maintain a slow spin rate. 
One idea is that the torques arising from the magentic 
connection between the star and disk can remove sub- 
stantial angular niomentum (e . g. iGhosh fc Lamb! 119781 : 
ICamenzind]lT990l: IKonigll Il99ll: Ishu et al.lll9'9l 7^When 



these torques are strong enough to enforce an equilib- 
rium stellar spi n rate, this idea is gen erally referred to as 
"disk locking" (jChoi fc Herbstlll996( ). The other idea to 
explain slowly rotating accretors is that powerful stellar 
winds are primarily resp onsible for removing angular mo- 
mentu m from the s tar (iHartmann fc MacGregoil 119821: 
;Mest elll984l: iHartma nn fc Stauffcr 19 89l:lTout fc PringL 



Paatz fc Camenzindiil996ciMatrfc Pudritdl200'5 



Numerical, dynamical simulations of the star- 
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disk intera ction (such as the investigations of spin equi- 
hbrium bv lRomanova et al]|2002l : [Long et aLl ;2005) typ- 
icahy exhibit significant torques arising both from the 
magnetic star-disk connection and from winds. 

The goal of the present paper is to further examine 
the possibihty of disk locking as an explanation for the 
slow rotation of young stars. To this end, we develop a 
model for the time-evolution of stellar spin rates, under 
the influence of magnetic star-disk interaction torques. 
In this work, we neglect the influence of stellar winds, as 
is customary in the disk-locking models. 

There are generally two types of di sk locking mod- 
els. One is the X-wind model (e.g.. iShu et al.l 119941 : 



lOstriker fc ShulfTOOl IMohantv fc Shull2008f) . This model 
is developed under the assumption that the star-disk sys- 
tem on average exists in a spin equilibrium state, in which 
the accreting star feels no torque. Thus, the X-wind 
model cannot be used to address whether or how the sys- 
tem reaches this equilibrium, nor how the system behaves 
when out of equilibrium (e.g., due to variability or due to 
the evolution of the star or dis k). All other disk l ocking 
models (e.g .,'Camcnzindlll993 IKonigllfToMt M^MM 
iMatt fc Pu dritz 2004) are of the second type, based on 
the general picture developed bv lGhosh fc LambI (|1978l ) 
for accreting neutron stars. The Ghosh & Lamb model 
provides a method for calculating the torque on the star, 
due to the star-disk interaction, for any state of the sys- 
tem. The star may spin up or down, depending on con- 
ditions, and there is a hypothetical disk-locked state, in 
which the net torque on the star is zero. Thus, the Ghosh 
fc Lamb model can be used to compute the evolution of 
accreting star spins and to assess how far from equilib- 
ri um a system may be. 

jCameron&Campbell" (TMS, hereafter CC93), |y1 
(|1994, 1991), and iArmitagc fc Clarke (1996, hereafter 
AC96) computed the evolution of pre-main-sequence star 
spins, adopting a Ghosh & Lamb-type model for torques 
on the star, and also including stellar contraction, a de- 
crease in accretion rate with time, and a prescription 
for the evolution of the magnetic field. For parameters 
that reasonably represent T Tauri stars, their models 
produced spin rates within the observed range. These 
models also showed that, regardless of the initial spin of 
the star, the spin rates rapidly (in less than ~ 1 Myr) 
approached an equilibrium value and stayed near equi- 
librium during the first few Myr of evolution. That is, 
the initial spin conditions were "erased," and the stellar 
spin was "locked" to a rate that depended only on the 
stellar mass and radius, magnetic field strength, and ac- 
cretion rate. These results generally supported the idea 
that disk locking explains the slow rotation of accreting 
T Tauri stars. 

However, since that work, it has been recognized 
that there is a serious theoretical problem with a key 
assumption of the classical Ghosh & Lamb torque 
model. The assumption is that the stellar mag- 
netic field lines connect to a large region of the disk 
and become highly twisted in the azimuthal direc- 
tion. This twisting is what gives rise to the mag- 
netic spin-down torque felt by the star. Several authors 
(e.g.. Ivan Ballegooiienlll 994l: iLvnden-BeU fc Boilv'11994!: 
Lovelace et al.lll995l: iHa vashi ct al. 1996; Miller & Stone 
1997!; 'Goodson et al." 1997; Fcndt & Elstner i2000t 
Matt ct al.,i2002, ; Uzdc nsky ct al., 2002; Romanova et al.l 



I2002t iKliker et al.1 120031: iBessolaz et all I2008D have now 
shown that, when the dipole field is twisted azimuthally 
past a threshold of approximately 45°, the magnetic pres- 
sure associated with the azimuthal component pushes the 
field outward. This leads to an infiation and opening of 
the field loops in a stellar rotation time scale, ultimately 
discon necting the star from the disk. lUzdenskv et al.l 
()2002|) developed a method for computing which mag- 
netic field lines would open. The amount of field opening 
depends primarily on how strongly coupled the magnetic 
field is to the disk, since the coupling is responsible for 
the twisting, and the twisting is responsible for the field 
opening. 

In or der to assess how the fi eld opening affects the 
torques. iMatt fc Pudritl ()2005U hereafter MP05) modi- 
fied the Ghosh fc Lamb formulation to in clude the effects 
of field opening, following the method of Uzdcnskv et aO 
(2002). MP05 showed that, for the relatively strong 
coupling expected in these systems, the opening of the 
field significantly reduces the spin-down torque on the 
star, relative to the classical assumption of a closed field. 
Thus, the equilibrium spin rate predicted by the disk 
locking picture is much faster, which calls into question 
whether disk locking can explain the existence of slowly 
rotating accretors. 

However, since the MP05 model could only calculate 
the torque for a given set of parameters (at a given 
epoch), their conclusion is based on the calculation of 
the equilibrium spin rate. MP05 were not able to de- 
scribe how the spin evolves in time. Therefore, there is a 
question as to whether the the stars actually evolve near 
equilibrium, and whether the evolution of the system, in- 
cluding changes in the accretion rate, stellar radius, and 
stellar spin rate may affect this conclusion. Also, it is 
important to develop a physical model for understand- 
ing the evolution of stellar spin over a range of ages. For 
these reasons, in this paper, we develop and utilize a 
stellar spin evolution model (similar to CC93; AC96:lYJl 
11994 I1995D that includes the updated torque theory of 
MP05. In this first effort, we do not attempt to explain 
all phenomena related to young star spins. Rather, we 
focus on testing whether models can produce spin rates 
within the observed typical range of ^ 1-10 days. This 
work extends the MP05 formulation to the time domain, 
during the first 3x 10^ yr (i.e., during the Hayashi phase). 

Section [2] contains a description of the model and de- 
tails of our calculations. Section [3] contains the results of 
stellar spin calculations that include the effect of mag- 
netic field opening. A discussion and summary of the 
conclusions of this work are contained in section 01 As 
a test of the model, and for illustrative purposes, we in- 
clude an Appendix containing the results of stellar spin 
calculations under the classical assumption of a com- 
pletely closed field. 

2. STELLAR SPIN EVOLUTION MODEL 

The goal of the present work is to determine the effect 
of different star-disk interaction torques on the evolu- 
tion of the spin rates of young accreting stars. For this 
purpose, in addition to calculating the torques, it is nec- 
essary to follow the evolution of the stellar mass (M*), 
radius (i?*), and moment of inertia, /» = k^M^R^, where 
fc^ is the normalized radius of gyration. This section de- 
scribes the complete model, assumptions, and adopted 
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Fig. 1. — Two different accretion histories considered in the mod- 
els, following equation |(TJ. Here and in all figures, a red line cor- 
responds to the case with the highest accretion rate considered, 
while a black line is for the lowest accretion rate. 

parameters. 

All of our models follow a one solar mass protostar 
as it evolves along the Hayashi (fully convective) track, 
before the formation of a radiative core. The "birthline" 
age of a one solar mass star is ~ 10^ yr (Stabler 1983) . 
which corresponds to the youngest ages at which stars 
are observed. To allow the system some time to respond 
to the initial conditions, we begin our computation at 
a somewhat earlier time of 3 x 10^ yr. We follow the 
evolution until an age of three million years, which is 
approximately th e end of the Hayashi track phase (e.g., 
ISiess et al\\2mh . 

2.1. Mass Accretion 

The torques acting on the star (detailed in i)2.4p arise 
from the interaction between the star and a surrounding 
accretion disk. As far as this interaction is concerned, 
one of the most important properties of the surrounding 
disk is the accretion rate. Specifically, this is the rate at 
which matter is fed into the interaction region (within 
a distance of several i?*), which we also assume is the 
rate at which material accretes onto the stellar surface. 
To follow the decrease of the accretion rate in ti me, we 
adop t an exponential decay (also used by CC93: IYiI|1994 

EMI), 

Md 



M„ 



-t/ta 



ta 



(1) 



where ta is the decay timescale and M^i is the "disk 
mass," equal to the total amount of mass that would 
accrete from time t — ^^ oo. We adopt ta — 10^ yr, for 
all models. 

Figure [1] shows the evolution of the accretion rate, 
given by equation ([1]). In all of our models, we con- 
sider two different disk masses, Md = 0.1 and 0.01 Mq, 
to sample a range of accretion rates. We chose these 
two disk masses to represent a relatively high and low 
accretion rate, shown as a red and black solid line in fig- 
ure [l] respectively. The true accretion histories of young 
stars are likely to be much more complex than equation 
(dl, but this is not yet well- understood. The advantage 
of our approach is that it is simple, and it samples the 
wide range of accretion rates that is observed in optically 
visible accreting young stars (e.g.. ,Hartmann et al.lll998l : 
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Fig. 2. — Evolution of stellar mass, normalized to solar units, for 
the two different accretion histories considered. The red and black 
lines represent the high and low accretion rates, respectively. For 
each case, we chose the initial stellar mass such that the final mass 
would equal one solar mass. 



'Sicilia-Aguilar et al.l[2005l iNatta et~alll2006[ ). 

Figure [2] shows the evolution of the stellar mass, which 
only depends on the prescribed accretion rate (eq. [1]) 
and an initial condition (at to = 3 x 10* yr). For each of 
the two different accretion rates considered here, we set 
the initial stellar mass such that the final mass equals 
one solar mass. That is, the initial stellar mass equals 
one solar mass minus Md, corresponding to 0.9 and 0.99 
Mq for the two cases. 

2.2. Stellar Structure and Evolution 

For the evolution of stellar r adius, we ad opt the sim- 
ple treatment of CC93 andlYil (|1994lll995D . This treat- 
ment models the structure of the star as a polytrope 
with an index of n = 3/2 and assumes a constant photo- 
spheric temperature, Te, during the Hayashi phase. The 
polytropic model results in a mean radius of gyration of 
fc^ = 0.2. The radiated, blackbody luminosity of the star 
is powered only by the release of gravitational potential 
energy, as the star contracts and accretes matter. The 
evolution of the stellar radius is then given by 



IT 



2— Ma 

M, 
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28TT<TRtT, 
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(2) 



where G is Newton's grativational constant and a is the 
Stefan-Boltzmann constant. 

The black and red solid lines in figure [3] show the evo- 
lution of stellar radius with time, resulting from equa- 
tions ([T]) and ([2]). In all of our models, we adopted an 
initial stellar radius of SRq and a photospheric temper- 
ature of Te = 4280 K. We chose these values in order to 
best match the radius evolution of a one solar mass star 
(with Z= 0.02) predicted by the more sophisticated stellar 
model of lSiess et al.l ()20Q0h . Also, the chosen value of T^ 
is within the narrow range of (nearly constant) Hayash i 
track temperatures exhibited by the ISiess et all (|20000 
mo del. For compari son, the circles in the figure show 
the ISiess et al.l (|2000f l model predictions. 

It is clear from figure [3] that there is a "bump" in the 
Siess et al. radius evolution from ~ 1-5 x 10^ yr that 
is not matched by our model. This is due to deuterium 
burning, which we neglect. Still, our simple treatment 
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t (xlO^yr) 

Fig. 3. — Evolution of stellar radius normalized to solar units. 
The solid lines show solutions of equation l[2}, with the black and 
red lines corresponding to the low and high accretion rate s (re- 
spectively) shown in figure [T] and discussed in section 12.11 For 
comparison, the c ircular symbols s how the prediction of a more 
detailed model bv lSiess et al.l II200CII '). 

matches the Siess et ah model at the begining and end of 
the evolution, demonstrating that the overall evolution 
of stellar radius during this phase is dominated by the 
Kelvin-Helmholz contraction. Thus, our simple treat- 
ment is sufficient for the present purpose. Also, the 
treatment of equation ^ captures the effects of mass 
accretion, which is neglecte d in more sophist icated mod- 
els (with the exception of iTout et al.lll99"9[ ). To com- 
pute the evolution for a longer time, beyond the Hayashi 
track, would require the use of a more sophisticated stel- 
lar model, which we leave for future work. 

2.3. Spin Evolution 

Along with the prescription above for the evolution 
of the star's mass, radius, and moment of inertia, the 
model also computes the angular momentum evolution. 
Assuming the star rotates as a solid body, the angular 
momentum equation can be expressed in terms of spin 
evolution, 













(3) 



where fi* is the angular velocity of the star, and T* is 
the net torque on the star, discussed in section [2^ 

It is sometimes convenient to express the spin rate as 
a fraction of the breakup speed, defined as the Keplerian 
velocity at the star's equator. This normalized spin rate 
is defined 



f = n. 






(4) 



In all of our models, we will consider two cases with dif- 
ferent initial spin rates. The two cases have initial frac- 
tional spin rates of /o = 0.3 and 0.06, representing the 
two extremes of rapid and slow initial rotation. In all of 
our models, we enforce a maximum spin rate of / = 1, 
since at this spin rate, some of the model assumptions 
(e.g., regarding the stellar structure) break down. As we 
show below, this limit is only reached in some of the most 
extreme models. 



2.4. Star-Disk Interaction Torques 

To calculate the torque on the star arising from ac- 
cretion and the magnetic interaction between the star 
and disk, we follow the formulation of MP05. Here, we 
summarize the parts of that work that are relevant for 
our spin evolution model, and the reader will find details 
in that paper. This torque m odel follows several pr evi- 
ous works ( e.g., CC93; AC 96; 'Ghos h fc La^[T97l |y1 
11994 119951: ILovelace et all [1995). but the primary ad- 
vantage the MP05 formulation has for this work is that 
it includes the effects of varying magnetic coupling and 
magnectic connectedness between the star and disk. 

2.4.1. Magnetic Field Prescription 

The basic model is one-dimensional (in radius, r) and 
assumes the star has a rotation-axis-aligned dipolar mag- 
netic field anchored into its surface, such that the mag- 
netic field strength varies in the equatorial plane as 



B, — _B* 



R. 



(5) 



where i?* is the magnetic field strength at the equator 
and surface of the star. The surrounding accretion disk 
is assumed to be thin, so that the magnetic field has a 
negligible radial component along the disk surface. 

The torque depends strongly upon the strength of the 
large-scale magnetic field i?*. Both disk locking and 
stellar wind spin-down models for T Tauri stars require 
global (dipolar) magnetic fields with surface strengths 
in the range of hundreds to thousands of Gauss. Ob- 
servations of pre-main-sequence stars typically find field 
strenghts with a mean absolute value of ^ 1-3 kG (e.g., 
Basri et al.l [1992: Johns-KruU 2007^ iJohns-KruU et alJ 



2009( ). while the large-scale (global) fields of these 
stars appear to be a t most several h undred Gauss 
JSaficr 1998; Johns-Kr Tni et al.l Il999al: iSmirnov et al 



lYang et al. 2007; Bouvier e t al.ll2007l: iDonati et al 



L2008; Hussain et al. 2009). While this discrepancy 
means that the magnetic fields are complex (not simple 
dipoles), it is clear that th ese stars possess dy namically 
significant magnetic fields (jGregorv et al.|[2008l ). For the 
models presented here, we choose two values of the mag- 
netic field strength, B* = 500 G and 2000 G, in addition 
to some cases with B^ — 0, used for comparison. 

Previous models for pre-rnain-se quence spin evolution 
(e.g., CC93; AG96: lYil ll99l 119951 ) chose a prescription 
for the magnetic field that depends upon the stellar spin 
rate and stellar radius, as motivated by the expectations 
of a stellar dynamo. By contrast, for the models pre- 
sented here, we choose to keep the magnetic field con- 
stant in tim e (a constant field was also conside red in 
the works of lJohns-Krull fc Gaffordir2002 and IyI [1991 
for three reasons. First, observations of mean fields of 
T Tauri stars find a relativ ely constant field strength 
(|Johns-Krull fc Gaff ord 20021: 1 Johns-K^[2007l) . regard- 
less of stellar age, radi us, or spin rate. Measu rements of 
large-scale field ( Johns-Krull et al. 1999a: Sm irnov et al. 



2004; Yang et al. 20 071: iBquvier et aljr20d7l : IDonati eTaL 



,2007. . , 2008;, Hussain et al.ll2009D are more rare, and there 
are not yet enough clear detections to look for trends. 
So a fixed field is consistent with observations. Sec- 
ond, the observed relationship between X-ray activity 
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and spin rates of T Tauri stars (jStassun at al.|[200^ sug- 
gests that they are pri marily in the "satu rated" or "su- 
persaturated" regime ([Johns- Krulll l2007f) , where solar- 
like stellar dynamo relationships are thought to break 
down. It is not yet clear how the global magnetic field 
behaves in this ful ly convective phase (Bro wning 2008t 
iDonati et al.ll20Q8l ). but at this point, adopting a con- 
stant field strength seems appropriate. Finally, the main 
purpose of the present work is to demonstrate how the 
stellar spin rate is influenced by the magnetic coupling 
to the disk. Thus, it is also convenient for illustrative 
purposes to keep the field a fixed constant. 

2.4.2. Magnetic Coupling to the Disk 

The stellar magnetic field lines vertically threading the 
disk are imperfectly coupled to it. To be clear, through- 
out this paper, the term "coupling" refers to the ex- 
tent to which the magnetic field is "frozen into" the disk 
material — that is, the effective magnetic diffusion in the 
disk. By contrast, the magnetic "connection" refers to 
the magnetic field loops that are attached both to the 
star and disk. MP05 showed how the strength of the 
coupling affects the amount of connected fiux, and how 
this in turn affects the integrated torques. 

The disk material rotates with Keplerian speed, which 
means that there is a singular location where the angu- 
lar rotation rate of the disk equals that of the star, the 
corotation radius, 



^rn = 



( GNU 



1/3 



= r^'^'R* 



(6) 



Every location, except for Rco, rotates at a different an- 
gular rate than the star. Thus, the magnetic field con- 
nected to the star and threading the disk will be twisted 
azimuthally (in the i/i-direction). It is assumed that the 
material above both the disk and star (the "corona") is 
of sufficiently low density (high Alfven speed) that the 
twisting of the magnetic field at the disk surface is equili- 
brated along the field line on a rapid timescale. Thus, the 
twisting of the magnetic field by the differential rotation 
between the star and disk imparts a torque, transferring 
angular momentum between the two. 

The variable 7(r) = B^/Bz quantifies the twist of 
the field at the surface of the disk at each radial loca- 
tion. This twisting occurs rapidly, on a stellar rotation 
timescale, so that a steady-state 7(r) is simply deter- 
mined by the balance between the differential rotation 
and the tendancy for the magnetic field to untwist by 
slipping (or reconnecting) azimuthally through the disk. 
This slipping rate depends on how well the magnetic field 
is coupled to the disk material. 

To describe the coupling, MP05 adopted a dimension- 
less magnetic diffusivity parameter, which they assume 
to be constant throughout the disk. 



hvk 



(7) 



where "qt is the effective magnetic diffusivity, h is the 
vertical thickness of the disk, and Vk is the Keplerian ro- 
tation velocity. In other words, (3 is the inverse of the 
effective magnetic Reynolds number of the star-disk in- 
teraction. A large value of (3 corresponds to weak mag- 
netic coupling (rapid field slipping through the disk) , and 



a small value to strong coupling (slow slipping). The 
various models of star-disk interaction torques in the lit- 
erature differ in their treatment of the magnetic field 
coupling, but all models adopt parameters that give sim- 
ilar results to MP05 with a value of (3 of order unitjQ. 
However, /3 is a highly unknown parameter, and MP05 
suggested that /3 < 10~^ (i.e., large magnetic Reynolds 
number) may be more reasonable for accreting pre-main- 
sequence stars. In our models presented here, we consider 
(3 = 10~^ to represent the most realistic case, but we also 
show cases with (3—1, for illustrative purposes. 

2.4.3. Magnetic Connection State and yc 

There is a dynamical limit to how strongly a dipo- 
lar magnetic field can be twisted a z imuthally. As 
show n by seve ral authors (lAlyl Il985l l Alv fc KuiipersI 



1990t iLvnden-B cll & Boilv 1994; Lovel ace et alT 



Bardou fc HevyaertSiil99& : iAe:apitou fc Papaloizoul 



1995 : 



20001: 



Uzdenskv et al.l [2002t ) . when the absolute value of the 
twist 7 exceeds a value near unity, the magnetic pressure 
force associated with i?^ pushes outward and leads to an 
opening of the dipole field loops. For practical purposes, 
this simply means that the star and disk are no longer 
causally linked, and no torques can be transmitted be- 
tween the two. This does not mean the open field lines 
necessarily impart zero torque. Open field lines can, and 
likely do, transport some angular momentum from the 
star and/or disk (via winds). However, no angular mo- 
mentum is exchanged between the star and disk along 
open field lines. In the present work, we consider only 
the torques arising in the star-disk interaction and thus 
neglect any torques that may arise from open field lines. 

To take into account the opening of the field, the MP05 
formulation includes a maximum for the absolute value 
of the twist, 7c. In the regions of the disk where the twist 
is greater than this critical value, the model assumes that 
the magnetic field no longer connects the star to the disk. 
The torque on the star that would have arisen from those 
magnetic field lines is instead taken to be zero. MP05 
showed that a value of 7c = oo represents the classical 
assumption of a field that remains connected at all radii 
in the disk. To show the effect of the opening of the 
magnetic field i n more realistic system s, a value of 7c = 1 
is appropriate (lUzdens kv et al.ll2002t ). 

Furthermore, MP05 showed that there exists a mode 
change in the magnetic connection between the star and 
disk, at a threshold value of the spin rate. Specifically, 
if 



/<(l-/37c)(7c^)-'/', 



where 



^ 



2B^R 



5/2 



MaVGM; 



(8) 



(9) 



then the stellar magnetic field only connects to a small 
region near the inner edge of the disk. In this case, which 
they call "State 1," no magnetic field connects the star to 
a radius in the disk that is larger than i?co- This is signif- 
icant, since in this state there exist no star-disk interac- 
tion torques that act to spin the star down (i.e., there are 

^ The one exception to this is the model of CC93 (see the Ap- 
pendix). 
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only spin- up torques). Alternatively, when equation ([S]) 
is not valid (for faster spin), the system is in "State 2," 
characterized by a magnetic connection to the disk over 
a range of radii on either side of Rco- In State 2, the star 
experiences both spin-up and spin-down torques. Fig- 
ure 3 of MP05 illustrates these two magnetic connection 
states 0. 

As evident in equation ([8]), magnetic connection State 
1 only occurs for relatively extreme cases of slow spin 
rate, weak field, and/or high accretion rate. As it turns 
out, all of the models with non-zero i?, presented in this 
paper remain in State 2 for the entire time simulated. 

2.4.4. Truncation Radius and Accretion Torque 

In the region of the disk at smaller radii than Rco, the 
twist of the field is such that torques remove angular 
momentum from the disk (giving it to the star). If the 
magnetic field is strong enough, there will be a radius 
above the stellar surface at which these torques can ex- 
tract all of the angular momentum of the material in the 
disk at that radius. Under these circumstances, the disk 
will be truncated, and accretion will proceed as a free-fall 
of material onto the star along the dipolar magnetic field 
lines. 

In State 1 [eq. (|S])], the location of the truncation radius 
is determined by 

(State 1) Rt - {-ic^f/^R*. (10) 

In State 2, the location of the truncation radius obeys 



(State 2) ( — 



-7/2 



1- 



3l 

Rco 



3/2' 



P 



i,pr 



m 



which we solve for Rt < Rco using a Newton-Raphson 
method. In our models, for each timestep, we first use 
equation ([5]) to determine the magnetic connection state, 
then either equation ([T0|) or ([TT|) to find Rt- Finally, we 
check that Rt is greater than i?,. If not, this means the 
magnetic field is not strong enough to truncate the disk 
above the stellar surface, and we reset Rt to equal i?*. 

The truncation of the disk is due to the stellar magnetic 
field extracting the angular momentum of disk material 
at a radius of Rt- Thus, the truncation of the disk and 
subsequent accretion of material adds angular momen- 
tum to the star at a rate equal to 



Ta = MaVGMJli, (12) 

which is referred to as the "accretion torque'fl Note 
that, although equation (fT2|) does not explicitly contain 
a dependence on the magnetic field strength, the mag- 
netic field is central to the determination of Rt (as long as 
Rt > R*)- Furthermore, when Rt > R^., material falling 

^ In the present work, we do not consider "State 3" 
discussed by MP05, whic h represents the "propeller" regime 
mUarionov fc !Junvaev|[l975l) in which there is no accretion of ma- 
terial onto the star. 

^ Equation (19) of MP05 for the accretion torque had an extra 
(negligible) term in brackets of k'^f, which is not correct. This 
term actually represents a term already in the basic angular mo- 
mentum equation, which is due to the change in the stellar moment 
of inertia (specifically, the 0«fc^R^Ma term). Equation II12I I of this 
work, which has also been derived b y many previous authors (e.g., 
CC93; AC96: [^hosh fc Lambiri979l) . represents the physical accre- 
tion torque. 



onto the star transfers the bulk of its angular momentum 
to the magnetic field before the material reaches the ste l- 
lar surface (e.g.. lRomanova et al.ll2002l : lLong et al.l20 05!). 
Thus, as long as Rt > i?*, the torque of equation ([T^ is 
experienced by the star as a magnetic torque, although 
we still refer to this as the "accretion torque" to distin- 
guish it from the component of the torque described in 
section 12.4.51 

2.4.5. Magnetic Torque 

In addition to the accretion torque, if the system is in 
State 2, the magnetic connection over a range of radii in 
the disk transports angular momentum between the star 
and the disk. In this case 



(State 2) T„i = 






[2(l + /37,)-i-(l + /?7,)- 



-2{Rco/Rtf/^ + iRco/Rt)% (13) 

is the net torque felt by the star from both spin-up and 
spin-down components. We refer to this as the "magnetic 
torque." If the system is in State 1, the magnetic field 
does not connect to a significant range of radii in the 
disk, so we set 



(State 1) T„ 



0. 



(14) 



As discussed in MP05, we assume that the disk will 
be capable of transporting away any angular momentum 
that it receives from the star (i.e., from negative torques 
in eq. [13]). For example, the disk may restructur e itself 
in re sponse to external torques (e.g., Rappaport et al.l 
|2004[ ) in such a way as to increase the effectiveness of 
viscous or turbulent angular momentum transport. 

2.4.6. Total Torque and Equilibrium Spin Rate 

The total torque on the star is the sum of the accretion 
and magnetic torques, 

n^Ta + r™, (15) 

and we neglect any other torques that may be present 
(e.g., from a stellar wind). The total torque can be either 
postive or negative, depending on the balance between 
spin-up torques (due to the trucation of the disk and 
to magnetic field lines that are azimuthally twisted in 
a direction that leads the rotation of the star) and spin- 
down torques (due to magnetic field lines that are twisted 
in a direction that lags the rotation of the star). 

For any given system, there is a particular stellar spin 
rate at which the net torque will equal zero (T* — 0). 
This is called the equilibrium spin rate, and the torque 
always acts to change the stellar spin in a direction to- 
ward this equilibrium. As indicated in equation ([3]), the 
spin also changes in response to changes in stellar ra- 
dius and mass, and in general to changes in the moment 
of inertia. If the contribution from the torque is large 
compared to that from changes in the moment of iner- 
tia, the stellar spin will evolve asymptotically toward the 
equilibrium rate. 

This general idea provides the theoretical underpin 



ning for the idea o f disk locking (iGhosh fc LambI 



Camenzindl 119901 : iKonigll Il99ll: lEdwards et a! 



1978 : 



19931: 



Shu et al.lll994[ ). which holds that the stellar spin rate 



is determined by the current conditions in the star-disk 
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Fig. 4. — Evolution of spin period (panel a) and spin rate expressed as a fraction of the breakup rate (panel b) for a star with no external 
torques (T* = 0). The four cases shown represent the two choices of accretion history and the two initial spin rates. The red and black 
lines correspond to the high and low accretion cases, respectively. The solid and dashed lines correspond to models with an initial spin 
rate of 0.06 and 0.3 times the breakup rate, respectively. 

interaction and independent of any initial conditions. 
MP05 described the method of computing the stellar spin 
rate in this hypothetical equilibrium state for their torque 
formulation, which uses equations (23)-(27) of that work. 
In presenting our results below, we compare the actual 
spin rates to the equilibrium values to gain insight into 
these systems. 

2.5. Numerical Method 

The coupled equations ([1]), ^, and ^ describe the 
evolution of the system. We wrote a computational code 
that solves these simul taneously, u s ing th e fourth-order 
Runge-Kutta scheme of lPress et al.l (|1994| ) , starting from 
ig = 3 X 10^ yr and ending at 3 Myr. For computa- 
tional efficiency and numerical stability, we use a dy- 
namic timestep. At each step, we compute the next 
timestep by requiring that none of the three main vari- 
ables, ri*, R^,, M*, change by more than 1% per step. 
This gives results converged to three significant figures 
and typically requires a few hundred timesteps for the 
entire evolution. 

In order to compute the torque at each timestep, the 
code first solves equation ([5]) to determine the magnetic 
connection state of the system. Then, depending on the 
state, the code solves either equation (fTU|) or (ITT]) to de- 
termine the disk truncation radius. Wc enforce a mini- 
mum value of Rt = R*. Finally, the code calculates the 
accretion torque using equation ((T^ and the magnetic 
torque using either equation (jTSl) or (fT4|). The torque is 
assumed to be constant during each timestep. We also 
enforce a maximum on the spin rate, corresponding to 
/ = 1 (see 



3. RESULTS 

This section contains results from several models. Ta- 
ble [T] lists the parameters for all cases, in order of their 
presentation and grouped by the figures in which the re- 
sults appear. All models have the same initial stellar 
radius (Si?©) and effective temperature (Tg = 4280 K; 
see fig. [3] and !j2.2p . Also, all models end with a mass 
of Mq, so that the initial mass depends on the accretion 
rate (see i j2.ip . For each case in table [1] we ran four mod- 
els, one for each combination of two different initial spin 



rates (see g2.3|) and two different mass accretion rates 
(see i j2.1[) . These parameters are intended to represent a 
range that is appropriate for accreting T Tauri stars. 

Section [3 . II includes two simplified cases with no mag- 
netic fields: one in which T* = 0, so that there is no 
star-disk interaction; and one in which i3* = so that 
there are no magnetic effects, just disk accretion. Section 
13.2! contains the main results, models that include mag- 
netic fields with realistic field line opening (7c = 1) and 
different values of /3, for two different field strengths. For 
comparison, and as a verification of our model and com- 
putations, the appendix contains results for cases with 
the classical assumption of a fully closed magnetic field 
(7c = 00) and P — 1. 

3.1. T,j = and i3* = Cases 

To begin, it is useful to examine a few simplified cases. 
The main goal is to understand individually how the stel- 
lar contraction and accretion alone affects the spin evo- 
lution, before adding the effects of magnetic fields. We 
will consider two cases here: one in which the net torque 
on the star is forced to be zero; and one in which there 
are no magnetic fields, so that the only torque on the 
star is due to the accretion of material from the disk. 

Figure m shows the evolution of stellar spin for the case 
of zero torque (T,^ = in cq. [3]). The figure shows both 
the spin period, P* = 2tt/VI^,, as well as the spin rate 
expressed as a fraction of breakup speed (eq. [1]), as a 
function of time. In this figure and in each subsequent 
figure, we show four different simulations, representing 
all combinations of two different accretion rates and two 
different initial spin rates. The red (black) colored lines 
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Fig. 5. — Results for B* = 0. Panel (a) and (b) show the same quantities as figure[4] The solid, dashed, red, and black lines in all panels 
have the same meaning as in figure |4] Panel (c) shows the net torque on the star for each case. In this panel, the dashed and solid lines 
of a given color lie on top of each other, since the spin rate does not affect the torque when B* = 0. The horizontal dotted line indicates 
T« = for clarity. Panel (d) shows the location of the truncation radius {Rt, thick lines) and the corotation radius (i?co, thin lines; eq. [5]), 
in units of the stellar radius. Here, Rt = R„ for all models, since there is no magnetic field to truncate the disk above the stellar surface. 



correspond to the highest (lowest) accretion rate consid- 
ered (see i)2.1l and figure [T]) . The solid (dashed) lines 
correspond to the slow (fast) initial spins considered (see 

ED. 

The spin evolution of the stars in figure |3] is due to 
contraction and the accretion of mass containing no an- 
gular momentum. After the initial time of 3 x 10^ yr, the 
cases with higher accretion rates (red lines) spin some- 
what more slowly than the cases with lower accretion 
rates. This is because the stars in the high accretion 
cases are gaining more mass, while conserving their an- 
gular momentum. Still, this effect is relatively minor, 
and it is clear that the spin evolution is dominated by 
the change in stellar radius. Specifically, the contrac- 
tion from the beginning to the end of the evolution (see 
figure [3|) , leads to a change in spin period by a factor 
of approximately 20, and a change in / by a factor of 
approximately 2. 

Next, we consider the case where angular momentum 
is added to the star from the accretion of disk mate- 
rial. Since the disk is assumed to be in Keplerian or- 
bit, it has high specific angular momentum relative to 
the star. In the simplified case where the star has no 
magnetic field, the disk material will extend all the way 
to the surface of the star, forming a boundary layer 
(|Lvnden-Bell &: Pring"Iill974[) . The torque on the star 



from the boundary layer is given by equation (|12p. with 
Rt = R* ■ This behavior is captured automatically in our 
model by setting B^, = 0. 

Figure [5] shows this case for the two different accretion 
rates and two different intiial spin rates. In addition to 
the evolution of stellar spin (panels a and b) , the figure 
shows the net torque on the star (panel c) and the lo- 
cation of the disk truncation radius and the corotation 
radius (panel d). Panel (c) indicates that the torque on 
the star is always positive (adding angular momentum to 
the star). The torque does not depend on the spin rate of 
the star but depends most strongly on the accretion rate 
(eq. [12]). Thick lines in panel (d) all overlap, indicating 
that the disk extends all the way to the stellar surface 
{Rt — R^,, a,s expected) in all of these cases. The thin 
lines in panel d show how the corotation radius (eq. [6]) 
changes with spin rate. 

It is clear from figure [5] how the accretion torque in- 
fluences the spin evolution of the star. For the lowest 
accretion rate considered, after 3 Myr, the star spins 
about 10% to 50% faster than in the T* = case (fig. 
[4|). But for the highest accretion rate considered, the 
stars spin much faster. The most extreme case, with a 
high accretion rate and fast initial spin (dashed red line) , 
reached the / = 1 limit enforced by our model (see H2.3P 
in approximately 1.5 Myr. In the case with the slowest 
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Fig. 6. — Results for B, = 500 G, 7e = 1, and /3 = 0.01, in the same format as figure[5] In panels (a) and (b), the dotted red and black 
lines show the equilibrium value ( i|2.4.6| l for the high and low accretion rates, respectively. In panel (d), the thin and thick lines of a given 
color and type are indistinguishable, since Rt is very close to iJco in all four cases. 



initial spin and lowest accretion rate, the spin period is 
1-2 days during the last few Myr. But the other three 
cases have spin periods shorter than one day during that 
time. In order to explain the observed spin periods of 
several days, for the range of accretion rates considered 
here, these stars must experience significant spin down 
torques. The next subsection presents cases that include 
the additional effects of torques that arise in the mag- 
netic star-disk interaction. 

3.2. Effect of the Magnetic Star-Disk Interaction 

The results in the previous section demonstrate how 
the spins of accreting stars are expected to evolve in the 
absence of magnetic effects. In this section, we include 
the magnetic torques that arise from the magnetic star- 
disk interaction. The torque calculation follows MP05 
(described in i j2.4p . and we adopt 7c = 1 in order to 
include the effect of the opening of the field that is ex- 
pected to arise due to the differential rotation between 
the star and disk. 

Figure [S] shows the results for cases with a dipole field 
strength of i?, = 500 G at the equator of the star, and 
(3 = 0.01. This field strength is similar to the strongest 
yet measured dipol e field on an accr eting young star, 
600 G for BP Tau (IDonati et al.|[200lFI . This value for 

* Note that IDonati et al.l 120081) report the polar field strength 



P was suggested by MP05 to best represent real systems 
(see Ji2.4.2p . The four models are shown in the same 
format as figure [5] and represent case 01 in table [T] A 
comparison of the upper panels of figures [5] and H] reveals 
that the magnetic star-disk interaction has very little in- 
fiuence on the spin evolution, in this case. In most mod- 
els, the presence of a 500 G dipole field results in spin 
rates that are 10-30% faster at an age of 3 Myr than the 
i?, = case. Only the model with low accretion rate and 
fast initial spin (black, dashed line) spins more slowly (by 
~ 10%) than the non-magnetic case. 

Panel (d) of figure [6] indicates that in all models, the 
disk is truncated very close to the corotation radius (the 
lines for Rt and i?co are indistinguishable). Thus, the 500 
G dipole is strong enough to truncate the disk, leading 
to magnetospheric accretion, in spite of the fact that the 
stellar spin evolution is not largely affected. 

The red and black dotted lines in panels (a) and (b) 
of figure [6] represent the hypothetical equilibrium (disk 
locked) spin rate (see i 32.4.6p for the two different accre- 
tion rates considered. It is clear that the three models 
that spin faster than the non-magnetic case are all spin- 
ning more slowly than their corresponding equilibrium 
rate. This means that the net torque on these stars is 
positive (acting to spin up the star). Furthermore, a 



of 1.2 kG, while S* represents the equatorial value (see ^2.4.11 1. For 
a dipole, B« is a factor of two lower than the polar field strength. 
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Fig. 7. — Results for B, = 2000 G, 7c = 1, and /3 = 0.01, in the same format as figure|6] 



comparison of panels (c) in figures [5] and [6] indicates that 
these three models have a larger net torque in the mag- 
netized (01) case. This explains why these models end 
up spinning slightly faster than the non-magnetized case. 
The model with low accretion rate and fast initial spin 
(black, dashed line) evolves near the equilibrium rate, 
resulting in a near zero torque in panel (c). However, 
this is simply due to a coincidence of the initial condi- 
tion with the equilibrium rate, since the torque is not 
strong enouth to maintain equilibrium within million- 
year timescales. 

In the disk-locking scenario, the magnetic star-disk in- 
teraction torques are strong enough to keep the spins 
near their equilibrium rate. The case with B^ — 500 G 
presented in the appendix (CI, figure[S]), where field line 
opening is ignored, does approximately follow the ex- 
pectations of disk locking. However, this is clearly not 
true for the 01 case (figure [B]), since the spins of all four 
models at all times depend upon the choice of initial spin 
rate. The main result of case 01 is that, when consider- 
ing the opening of the field and strong magnetic coupling, 
a dipole magnetic field strength of 500 G is not strong 
enough to significantly infiuence the spin evolution of one 
solar mass accreting pre-main-sequence stars. 

Figure [7] shows the results for the 02 case, which is 
the same as 01, except that the magnetic field is much 
stronger {B^, — 2000 G). In the 02 case, the magnetic 
star-disk interaction has more of an influence on the spin 



evolution. By comparing the top panels of figures |6] and 
[71 it is clear that all models are spinning more slowly in 
the 02 case. 

The spin rates of the two models with the high ac- 
cretion rate (red, solid and red, dashed lines) converge 
toward the equilibrium value at an age of ~ 1 Myr. Thus, 
these models approach the disk locked state, and the ini- 
tial condition is "erased" at that time. The same is not 
true for the models with the low accretion rate (black, 
solid and black, dashed lines), since the spin rate after 
3 Myr still depends upon the intial spin rate for those 
models. 

The 02 case is significantly different than the case with 
i?* = 2000 G presented in the appendix (C2, figure ITO)) . 
which ignores field line opening. By comparison, the disk 
locking in the 02, high-accretion models is marginal and 
only happens after ~ 1 Myr, which is approaching the 
lifetime of disks. Furthermore, it is important to note 
that the rotation periods of all models in the 02 case are 
less than ~ 3 days in the age range of 1-3 Myr. This is 
consistent only with the fastest rotators in the observed 
spin distributions (see Sjl]). 

Thus far, we have considered /3 = 0.01, since this was 
suggested by MP05. However, MP05 showed that the 
spin-down torque felt by the star is the strongest for a 
magnetic coupling strength of /3 = 1. This is indeed a 
highly diffusive situation, as this means that a magnetic 
field line that is tilted by 45° will slip through the disk 
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Fig. 8. — Results for B, = 2000 G, 7c = 1, and /3 = 1, in the same format as figure|6] 



at a speed equal to the Keplerian speed (see MP05). Al- 
though this is unexpectedly diffusive for an astrophysical 
plasma, the correct value of (3 for these systems is highly 
uncertain. Thus, it is instructive to look at this case, 
which will have the most extreme effect on the stellar 
spins. 

Figure \8\ shows the results of this extreme case, with 
/3 = 1 and B^ = 2000 G (case 03). It is clear that, in 
this case, the star-disk interaction completely dominates 
the spin evolution of the stars. The spin rates converge 
toward their respective equilibrium values in 1-3 x 10^ yr. 
Subsequent evolution closely follows the equilibria, con- 
sistent with the disk locking scenario. The spin periods 
are in the range of 3-10 days for ages of 1-3 Myr. It is 
evident from panel (d) that the disk truncation radii are 
close to the corotation radii, but not so close that the 
lines overlap. This difference from the previous 2 cases 
is due primarily to the much higher magnetic diffusion 
rate in the 03 case. 

4. DISCUSSION AND CONCLUSIONS 

The star-disk interaction model of MP05 includes the 
effects of the field opening expected to arise from star- 
disk differential rotation (7c — 1), as well as the effects 
of magnetic coupling, parameterized by /3. Based essen- 
tially upon the idea that the disk likely exists in a state 
of high magnetic Reynolds number, MP05 argued that 
a value of /3 = 0.01 is appropriate for real systems. In 



the present work, we have extended the MP05 torque 
model into the time domain, in order to explore the con- 
sequences of magnetic star-disk coupling and connected- 
ness on the evolution of stellar spins. 

To this end, we developed a stellar spin evolution 
model (31]) • The model follows an accreting, one solar 
mass star during the Hayashi phase (from 3 x 10'' yr to 
3 Myr) and considers the torques expected to arise from 
the magnetic star-disk interaction. The primary goals of 
this work were to examine the role of disk locking and 
to determine if the star-disk interaction torques alone 
are sufficient to explain the observed range of spin rates. 
In developing our torque and spin-evolution model, we 
have adopted many of the same assumptions as the disk- 
locking models in the literature, except that the MP05 
torque formulation includes the effect of the opening of 
the magnetic field. Thus, while there are inherant un- 
certainties associated with the adoption of many of these 
assumptions, our results serve best to highlight the ef- 
fect of field opening relative to previous results in the 
literature. 

Given the expectation that 7c = 1 and (3 = 0.01 best- 
represents the conditions in real systems, and for a range 
of accretion rates and field strengths appropriate for T 
Tauri stars, the two main conclusions of the present work 
are as follows: 

1. The models presented in figures [6] and [7] exhibit spin 
periods ranging from 3 days to less than 1/3 day, in the 
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age range of 1-3 Myr (see Wd.2^ . These are consistent 
only with the fastest rotators in the observed spin pe- 
riod distribution, which is significantly populated from 
approximately 1-10 days (see [JT]). It is apparent that 
the torque arising from the magnetic connection between 
the star and disk is not sufficient to explain the relatively 
slowly spinning, young stars. 

2. Furthermore, the stars in these models are gener- 
ally not in spin equilibrium (with a net zero torque) , and 
the spin rate at all times depends upon the initial spin 
rate (see t j3.2[) . Only the models with the strongest field 
strength and highest accretion rate neared their equilib- 
rium spin rate in ~ 1 Myr (see figure [7]). However, this 
equilibrium had a relatively short spin period of less than 
a day. It is apparent that disk locking does not play a 
strong role, except possibly for the most rapid rotators. 

Note that the term "disk locking" generally refers to 
the idea that the angular spin rate of the star is nearly 
equal to that of the inner edge of the disk (i.e., Rt « 
Rco)- In our magnetic models, this condition was always 
true. In fact, the truncation radius was nearest to the 
corotation radius in the models that were furthest from 
spin equilibrium (e.g., compare panel (d) in figures[6]and 
[9]). Thus, these stars do have an angular spin rate that 
is very close to that of the disk inner edge, although 
it is not appropriate to think of these systems as being 
"locked" to any particular spin rate. Instead, this is a 
consequence of how the disk is truncated. When the 
magnetic coupling in the disk is strong (small /3) , the disk 
will generally be truncated close to the corotation radius 
(in State 2; see ^2.4.4p . But, the condition that Rt ~ i?co 
does not necessarily mean that the star experiences a net 
zero torque. 

In order to explain the existence of slowly rotating 
young stars, it is necessary to consider other possibili- 
ties. Since the magnetic torques depend strongly on the 
magnetic field strength, it is tempting to suggest that 
a stronger magnetic field may solve the problem. Us- 
ing our model, we find that (for /? = 0.01) a dipole field 
strength of _B* = 10'* G is required to maintain spin pe- 
riods consistent with the slow rotators. However, no T 
Tauri star has been found to have surfa ce field strengths 



greater than 3 kG () Johns- Krulll 120071 ). an d the global 
(dipole) fields are ge n erally even weaker ( e .g..lSafi 
iBouvier et all [20071 iDonati et all [20071 120081 ). Given 



d the global 
JSafieilll998l : 



the relatively small number of stars for which there are 
magnetic field measurements, the present situation could 
be improved by more and improved observations of the 

global field strength and geometry. 

Since the models presented in section 13.21 with weak 
magnetic coupling (figure [8]) exhibit slow spins, it is also 
appropriate to consider the possibility that real systems 
have weak coupling. As discussed above, the conditions 
that are expected to best represent young star-disk sys- 
tems from "first principles" suggest values of /3 = 0.01. 
However, /? is a highly uncertain parameter because the 
physics of magnetic coupling is not well understood, nor 
are the conditions that influence the coupling (e.g., the 
ionization fraction of disk gas, turbulence levels in the 
disk, or magnetic reconnection rates). In order for mag- 
netic coupling to explain the slow rotators, the coupling 
strength must not be very different (neither larger nor 
smaller) than (3 — 1, since MP05 showed that the max- 



imum spin-down torque occurs for a /3 of unity. The 
condition that /3 « 1 seems to require some fine-tuning, 
which makes this possiblity even more difficult to justify 
at present. 

Therefore, it appears necessary to consider effects that 
are alternative or additional to the magnetic torques 
arising from the star-disk connection. Real systems seem 
to be variable on all timescales (e.g.. [Hartmannlll997fl . 
and we have neglected variability here. It may be pos- 
sible, for example, that short timescale (e.g., < lO'* yr), 
large-magnitude variations of the accretion rate are 
important for t he stellar mass and angular momentum 
evolution (e.g.. [PophamI [l996[ ) . Finally, it is necessary 
to consider the angular momentum loss from pre-main- 
sequence stellar winds (see flj- The idea that stellar 
winds may be import ant during the accretion phase is 



well-s up ported (e.g.. iDecampli 119811 : 



weii-s up pc . .-_ 

1983 i Kwan fc Tadem ard [igSSf 
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Hartmann et al 
1,990; Fc ndt et al.l 1199 5: Fcndt & Camcnzind 199 
Safier, il998t Bogovaloy fc Tsinganoa i2001i: iSauty et aTl 



200a Edwards et al.l 120031 120061: iDupree et all 120051: 
Meliani ct al.' '2006'; Kurosawa et al.l l2006t iKwan et al.l 
2007; Fcndt 2009). Stellar wi nds may be powered 
by the accret ion process itself ([Tout fc Pringld 119921 : 
lCranmeill20 08t) and be the key driver of angular momen 
tum losspartmann & MacGrcgor 1982; Mcstcl 1984; 



Hartmann & StaufFcr 1989; Pa atz fc Camenzindi [1996; 
Matt & Pudritz 2005a, 2008a,b). 

Some recent magnetohydrodynamic (MHD) s imula- 
tions of ra pidly rotating stars (e.g., iRomanova et al.l 
120051 12009( ) exhibit a type of "propeller" regime in which 
intermittent accretion occurs on typical timescales of sev- 
eral orbits of the inner disk, while most of the would-be 
accreting material is instead launched into a wind. For 
the calculations in the present work, the torque repre- 
sents that which is averaged over a calculation timestep. 
Our timesteps range from approximately 400 years (at 
the earliest times) to ^ 10^ yrs. While some of the 
simulations in the work cited above have been run for 
thousands of dynamical times, this is still orders of mag- 
nitude shorter than the typical timestep in the present 
work. It is not clear whether the torque formulation we 
adopted accurately describes the time-average behavior 
of the accretion and closed field region in the propeller 
regime simulations. In addition, while these MHD sim- 
ulations reliably demonstrate many of the complexities 
and sometimes episodic nature of the magnetic star-disk 
interaction, when considering the long-term evolution of 
the star and disk, it is not yet clear whether the dura- 
tion of the episodic (e.g., propeller) regime is significant 
for the overall spin evolution of young stars. However, 
the MHD simulations such as those cited above often ex- 
hibit significant outflows and torques arising from open 
magnetic field regions. In this respect, our result that 
the torques arising only from closed field regions are not 
sufficient to significantly spin-down the star, agrees with 
the simulations. 

It is clear that the opening of field lines is an impor- 
tant effect, since including this in the models indicates 
considerably different physics at work than the closed 
field models, when comparing to the same observational 
data. In future work, we will extend the present analysis 
to include the angular momentum loss from accretion- 
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Fig. 9.— Results for B, = 500 G, -jc ■ 



powered stellar winds. 



and /3 = 1, in the same format as figure |6] 
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APPENDIX 
CASES WITH CLOSED MAGNETIC FIELD (7^ = 00) 

In order to demonstrate that our model and code reproduces the results of the classical models, and for comparison 
to the cases presented in section |3l here we present cases with 7c = 00. With the MP05 torque formulation, this 
value of 7c keeps the magnetic f ield connected e verywhere to the disk, regardless of how twisted the field becomes. 
This well- represents the classical IGhosh fc LambI (|l978) model, and results in stronger spin-down torques on the star. 
Most star-disk interaction and disk-locking models follow the Ghosh & Lamb model and assume a magnetic coupling 
strength that is equivalent (within ~20%) to a value of /? = 1. Thus, we adopt this value of (3 here. 

In table[Tl the two cases presented here are labeled CI and C2, which have i?, = 500 G and i?* = 2000 G, respectively. 
Figures |9] and \TU\ show results from these two cases, where the line styles and colors have the same meaning as in 
figure |6l The spin evolution of all models in both figures is dominated by the magnetic interaction between the star 
and disk. For the B* = 500 G case (figure [9]), the spin rates converge toward their equilibrium values in 2 x 10^ and 
10^ yr for the cases with high and low accretion rates, respectively. For the B* = 2000 G case (figure [10]), convergence 
occurs for all models within ^ 10^ yr. Thus, when the magnetic field is not allowed to open as expected, the evolution 
follows closely with the expectations of the disk locking scenario. 

As in all other magnetic models presented in this paper, the disk truncation radii remain near, although slightly 
inside, the corotation radii. However, in these cases with jc — 00, Rt is generally further away from Rco than in the 
cases presented in section [3] This is a consequence of the stronger spin-down torques and higher value of /3 for the 
models of this section. 

Examination of figures |9] and [10] reveals that the spin rates are slower when the magnetic field strength is higher or 
the accretion rate is lower. Under the conditions of disk locking, the angular spin rate of the star is proportional to 
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Fig. 10. — Results for B, = 2000 G, ■jc = oo, and /3 = 1, in the same format as figure|6] 

Ma S* , which is exhibited in the models presented here. Furthermore, the spin periods of ah the models are in the 
range of 1.5-20 days during the ages of 1-3 Myr. These values are consistent with the range of observed spin periods 
(see ^, and this has been the primary physical justification for the disk- locking picture being applied to young stars. 

The mode ls presented in this section produce very similar quantitative results to the spin evolution models of ml 
(|1994 I1995D and AC96. The main differences between all of these works is the treatment of the magnetic field and 
accretion rates. In particular, other works considered a magnetic field strength that evolved with stellar spin rate and 
stellar radius. Here, we have assumed a constant surface magnetic field strength. Our assumption results in different 
behavior from those models, particularly at times earlier than ~ 1 Myr. At these younger ages, the stars have relatively 
large radii, and the assumption of a constant surface field strength means that the total magnetic energy in the field 
(proportional to B^R^) is larger. This is the m ain reason that the equilibrium spin rates are slowest at earlier times, 
in our models. The models of |y3 ()1994 119951 ) and AC96 exhibit spin periods that are more ne arly constant in time, 
due the different prescription of the field. Still, when our results and those ofJY^ (J1994I . [19951 ) and AC96 are scaled 
according to the torque theory to similar values of accretion rates and field strengths, the results are quantitatively 
consistent. 

The spin evolution model of CC93 is somewhat different fro m other wo rks. The assumptions in that work lead to 
a magnetic coupling that is significantly stronger than in|Y3 ()1994 I1995D . AC96, and most other models. We have 
determined that we get similar results to CC93 for models with 7^ = c» and a coupling strength of /3 = 0.05. This 
relatively strong coupling leads to a large azimuthal twisting of the magnetic field. This, together with the requirement 
that the magnetic field remain closed, leads to very strong spin-down torques on the star. For this reason, the CC93 
model pred icts slower equilibrium spins than other models for a given field strength (e.g., see iJohns-Krull et al.lll999bl : 
lJohns-Krul]| |2007D . However, as shown in section |31 when the opening of the field is properly taken into account, the 
spin-down torques are weaker, and the magnitude of this effect is larger for cases with strong magnetic coupling (as 
in the CC93 model). 
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